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ABSTRACT 


In the magnetosphere, the equation of heat flux parallel to the magnetic 
field must be modified because of the low electron density and large mean free 
path. This leads to a density-dependent heat conductivity which is decreased 
such that the high temperatures in the magnetosphere observed by Serbu and 
Maier can be explained. In the upper magnetosphere, turbulent heat transfer 
in accordance with Bohm's diffusion coefficient can account for heat fluxes on 
the order of 10 erg/ cm sec across the magnetic field. Such heat fluxes are 

consistent with ion-electron temperature differences of Tj - T e ~ 10 4o K also 
discovered by Serbu and Maier . This direct thermal energy input is negligible 
for the energy balance of the magnetosphere when compared with the heat fluxes 
of 10 - 10 erg/cm sec generated through non-local heating by fast iono- 

spheric electrons. 
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THE MAGNETOSPHERIC TEMPERATURE DISTRIBUTION 

by t 

H. G. Mayr* and H. Volland 1 

Goddard Space Flight Center 


INTRODUCTION 


Measurements of the low-energy electron spectrum by Serbu and Maier 
[1966] on IMP II have shown that in the magnetosphere these electrons have an 
almost Maxwellian velocity distribution. 

In the equatorial plane the electron temperature T e increases with dis- 
tance r from the earth's center as 

T e °c r 2 , ( 1 ) 

reaching values of the order T e ~ 20, 000°K at heights between 5a and 15a where 
a is the earth radius. Otherwise relatively low temperatures of the order of 
3000°K have been observed at the base of the protonosphere at heights of about 
1000 km above the ground at low and high latitudes [Brace, Reddy and 
Mayr . 1967]. 

According to the classical formulas of heat conduction within an aniso- 
tropic plasma, the predominant heat flow parallel to the geomagnetic field 
within the protonosphere is due to electrons [ see, e. g. , Kaufman , 1966]. The 
coefficient of heat conductivity of electrons is 

k „ ~ k = AjT e 5 /2 erg/cm sec deg, (2) 
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with 


A, ~ 1 X10' 6 . 

Perpendicular to the magnetic field the ionic component predominates and 
its coefficient of heat conductivity is 


K 






a 2 


Nj 2 

B 2 T^/2 


erg/cm sec deg, (3) 


with 

A 2 ~ 3 X10' 16 . 

Here Tj, Nj, and are temperature, density, and collision number of 
the protons; B is the earth's magnetic field, the ion gyrofrequency, m^ the 
ion mass, and m e is the electron mass. 

Taking mean numerical values of T e ~ ~ 10 4 °K, ~ 10 2 cm -3 and 

B ~ 10 -3 gauss, leads to a ratio between the coefficients of heat conductivity of 

- = 3 x 1(T 12 ; (4) 

K „ 

thus, one must expect an extremely large anisotropy of temperature: namely, 
near-isothermy parallel to the field lines and almost complete heat-insulation 
perpendicular to the field lines. 


To derive a first approximation of the temperature profile within the pro- 
tonosphere, we consider a heat-insulated field line and neglect the geometries 
of the dipole field. The static energy equation has then the form 

div j = - div («grad T) = Q, (5) 

whereT is the heat flux vector and Q is the external energy input per unit volume. 
Introducing the coefficient k u of (2) into (5) and taking Q = 0 leads to the solution 
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T = 


7/ 7 k T^T 

t i /! + \ T7 


(S - St) 


2 /t 


C s 


Z A , 


( 6 ) 


where Tj and Tj are the lower boundary values at height Sp From Equation 

6 follows the expected small temperature increase (as s 2//? ) parallel to the field at 
greater heights where (s » Sj). To maintain even this small gradient and to reach 
a temperature Tj = 20, 000 e K at the top of the field line, a heat flux of the order of 

jj = KjTj = l erg/ cm 2 sec 

had to flow into the thermosphere at medium and higher latitudes. This would 
imply a temperature gradient of the order 

T^ = 100 e /km 

at a height of 1000 km. Tamao [1966], taking into account the geometry of the 
field lines, arrived at a similar result. He assumed that this heat flow must 
be supplied from interplanetary space via turbulent heat transfer across the 
field lines at greater heights. But the observations made by Brace, Reddy and 
Mayr [1967] suggest an upper limit of the electron temperature gradient at 
1000 km that is smaller by at least two orders of magnitude: 

Tj ~ l°/km. 

If we introduce this value into (6) together with Tj = 3000°K, we obtain a tem- 
perature at the top of the field lines of T II ~ 7000°K at a distance r = 8a. This 
temperature is considerably smaller than that obtained by Serbu and 
Maier [1966]. 


The following investigation seeks to resolve this striking discrepancy. A 
more rigorous treatment of heat conduction shows that the heat conductivity is 
effectively decreased in regions of low electron densities, causing a general 
increase of temperature at greater heights. 

Furthermore, the effect of turbulent heat transfer perpendicular to the 
magnetic field lines is estimated. This leads to an understanding of the 
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ion-electron temperature difference observed by Serbu and Maier [1967] in the 
upper magnetosphere. 


PARALLEL HEAT CONDUCTION 

A characteristic quantity for heat conduction is the mean free path. This 
quantity is given for electrons by 

T ^ 

X = K ~ 3 x 10 4 (7) 

and is ordinarily related to the parallel heat conductivity through 

(e) _ 

“ cXNg, (8) 

where c is the mean thermal velocity (c oc T e ^ z ). 

Equation (8) is valid only if the temperature variation can be considered 
linear within the range of the mean free path, and if the mean free path is 
almost constant over the same range; only then does (7) have real significance. 

In the magnetosphere (where N can be on the order of 10 1 - 10 2 cm~ 3 and 
T e can be on the order of 10 4o K) the mean free path (as defined in (7) ) is on the 
order of 10 4 - 10 6 km. (This distance is comparable with or exceeds the 
dimensions of the plasma confined in field tubes. ) In the light of Serbu and 
Maier's measurements it is obvious that the temperature variation is no longer 
linear over such a distance and that the mean free path is strongly variable. 
Thus, (8) is no longer valid, and we must modify the heat-conduction equation 
so as to take into account this new situation in the magnetosphere. 

In this paper, we approach the problem by means of a modified mean- 
free-path method. This is not precisely a rigorous treatment; rather, it seeks 
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to prove that heat conductivity becomes density-dependent and is marked by a 
decrease in the magnetosphere. 


Expanding the temperature into a power series in s around s 0 = 0 gives 

T = T 0 + (fr) o S + 2" (fif) o s2 + • • (9) 

where higher terms are assumed to be negligibly small. If £ is the angle of the 
velocity direction toward the s axis, each electron carries the energy 
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~2 kT o - *- cos i 
+ fer 



( 10 ) 


through the plane s = 0. The quantity is considered the mean free path of 
the electrons in the region where s<0; k is Boltzmann’s constant. Corre- 
spondingly, electrons going in the opposite direction transport the energy 


-§- kT Q + \ + cos f 
+ ¥- cos! e 



• > 


(ii) 


where \ + is the mean free path in the region where s > 0. The difference of 
these energies is 
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( 12 ) 


- ^\+ + X-) cos £ 


Nil 


- ~ (X+ 2 ~ X- ) COS £ 


[fe (i kT )] 


During unit time c cos £ N e fd £ electrons traverse the unit area at 
s = 0, if f is the distribution function for the velocities c, given by 


" 25 (wfne E T Y exp (Sr) 


m e 3 c 2 sin £ dc 


The net energy transported by electrons with velocity c and angle £ is therefore 


-(X+ + X-) c |j kT^ Netcos 1 tdt 




(i kT ) Ni 


(IS) 


e fcos 3 £ d £ 


Integration over 0s£s| and 0 s c s«o leads in a straightforward manner to heat 
conduction flux in the form 


j =4c (X ++ x.) N e |I-||c (4-XJ) N e f£. 


8s z 


(14) 


If we assume that the mean free path is almost constant in space, then A + -X. , 
and (14) reduces to the classical formula 


= -^f- XNp |l =-at 5 / z |X 


9s 


as 


(15) 


where X is defined by (7). This shows that the additional term in the expansion 
of the heat-flux formula (equation 14) arises from the variability of the mean 
free path. 
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The definitions of X + and X_ are rather problematic. A reasonable defi- 
nition of these quantities can be given in the form 


f jk „ f°ds_ _ 

X J-\. X 


where \ is defined by (7). This leads to 


M 


kf-r , ds - i/ N k* (g) • + J (0) -p 

* / o L 7 o v 'o J 


ds - 1, (16) 


and 


ds = 1 • (17) 


•A- -A. J 

A treatment of the mean free paths in this form would lead to the functions 


= X. 


Ar I!. . 9 4 T\ 

’ 8s W ) 


(16a) 


and 


. y /rr dT , a*T\ 

x - = x - w) 


(17a) 


Introducing the heat flux, (14), into the energy continuity, (5), and integrating 
up to the top of the field line (equator), leads to the differential equation 


I c (X + + X-)Ne |f + fj[ 6 Ne0 - 


(18) 


3 

- 7 'n + x -> N e„ (|f) + i a n' x +- xJ >n % fe) + / Q ds - 

x 'll x 'n s 
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Because of symmetry with respect to the equatorial plane, at s — s^ T at the equator 
the conditions 


= 0 

Wn 


(19) 


and 


x +n - x -n 


must be fulfilled. This yields the equation 

2 /* s n 

^ c (X + + X_) N e + “§3T ® ^ x + 2 “ x - 2 ) N e = J Q ds > ( 20 ) 

s 

for which a solution can be obtained if the temperature is known at some loca- 
tion Sj (T = Tj). 

Thus the problem is uniquely defined. The expansion of the heat flux in a 
form, (14), though it involves also second temperature derivatives, does not 
require an additional boundary condition for the solution of the energy equation. 


A solution as outlined here is difficult, because of the complicated func- 
tions X + and A_ (equations 16a and 17a). Therefore, we shall discuss a very 
specific magnetospheric condition that allows a rather simple approach. 

In the heat-flux formula, (14), the second term, proportional to the second 
temperature derivative, is negative. This is true because the magnetospheric 
temperature distribution exhibits a maximum at the equator; and because X + is 
larger than A_ for the observed positive temperature gradients in the magneto- 
sphere. That term in (14) therefore decreases the heat flux and thus decreases 
what might be called the "effective conductivity." For a constant temperature 
at the base of the magnetosphere, this decrease of the heat conductivity leads 
to an increase of the electron temperature in the upper magnetosphere. For 


8 


the sake of simplicity, we shall disregard this second term in (14) in our sub- 
sequent discussion; therefore, the derived temperature should be regarded as 
a lower limit. 

From (16) and (17) it is evident that at the bottom of the magnetosphere, 
in a region where N ~ 10 2 cm 3 and T ~ 4 X 10 3 °K 

^ + >5 x 10* km, 

4 

A - < 5 x 10 km, 

and consequently 

(*- + + ^_) > 5 x 10 km. 

Distances of this order are comparable with the length of field lines which, for 

4 

L = 6 , are about 5 x 10 km from the earth surface up to the equatorial plane. 
Because of this, the upper limit of (X + + \_) cannot exceed the total length of 
the field line — which is 10 5 km; it follows that 

10 5 > (X + +XJ> 5 x 10 4 km. (21) 

As can easily be verified, this condition is also valid along the entire field line 
above the base of the magnetosphere. 

Thus, it seems that under magnetospheric conditions of low densities the 
quantity \ + + is almost constant. Its lower limit is close to the dimensions 
of the field tubes and its upper limit cannot exceed this. The heat flux in the 
magnetosphere then takes the form 

j kc f N e 9T_ — Ml N 6 t' 4 il (22) 

4 3s 3s 

where, after substitution for c, M is a constant and S. is the length of the 
field line. 

A comparison of this modified heat-flux formula with the classical ex- 
pression (equation 15) shows marked differences. The classical formula gives 
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the conductivity as independent of density and increasing strongly with increas- 
ing temperatures. Equation 22 states that under magnetospheric conditions 
the heat conductivity becomes density-dependent (decreasing with density), 

and that the temperature dependence of the conductivity becomes very weak 
l/ 2 54 

(ccT ) as compared with the classical conductivity (ocT ). Thus, in the 

magnetosphere the conductivity, because of its density-dependence, tends to 

decrease with height; in striking contrast to the behavior of the classical 

conductivity which, because of its temperature dependence, strongly tends to 

increase with height. 


Combining the heat-flux equation 23 with the energy equation 5 leads to 


which yields 


MfNT 1 / 2 9_T 
as 


T 3 /2 = T 3 /2 + 3 1 

1 2 Mi 



(23) 


^ -12 

If we assume that T^ = 4 X 10 °K and Q~3 X 10 — a heating rate expected 

from fast electrons that escape the ionosphere [ Geisler and Bowhill , 19651 -and 
that l =10 , then we obtain a temperature distribution as shown in Figure 1. 
This reveals a magnetospheric temperature structure that, at greater heights, 
is consistent with Serbu and Maier's measurements of temperatures on the 

4 

order of 2 X 10 °K. Figure 1 is also consistent with the expected temperature- 

gradients at lower altitudes (on the order of l°K/km, corresponding to energy 

-2 2 

fluxes on the order of 10 erg/cm sec). 


We have shown that under magnetospheric conditions of low density the 
heat conductivity becomes dependent on density and almost independent of 
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FIELD LINE LENGTHS (KM) 



Figure 1. Calculated temperature distribution 

along a field line in the magnetosphere. 
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temperature. As a consequence the heat conductivity is significantly decreased; 
and from this the observed unexpectedly high temperatures in the magnetosphere 
can be verified. The increase in equatorial temperature could thus be explained 
by the pronounced decrease of electron density observed by Carpenter and 
Smith [1964], Serbu and Maier [1966], and Taylor, Brinton, and Smith [1965]. 

TURBULENT PERPENDICULAR HEAT TRANSFER 

D. Bohm [1949] found, first from experiment and then from theory, that 

the effective perpendicular diffusion coefficient is enormously enhanced in a 

turbulent plasma. This coefficient D is independent of density and related to 

B 

temperature T and magnetic field B by the formula 

Dg ~ 5.4 x 10 2 T/B r Spitzer 1956] (24) 

The coefficient of turbulent heat transfer is proportional to D : 

B 

Kg=kND B . (25) 

At magnetospheric heights of 10 earth radii, this value is almost 10 7 
as large as the classical coefficient of heat transport based on collisions 
(equation 3). Combining (3), (24), and (25) (and equating T. to T and N. to N), 
gives 



BT 3 / 2 

N 


(26) 


Assuming Bar 3 and Serbu and Maier 1 s [1966] experimentally derived tempera- 
ture and density variations, 


T oc r 2 , N oc r" 


(27) 


12 


gives 



(28) 


As we go farther from the earth, the heat-transfer coefficient based on 
Bohm's diffusion coefficient increases strongly as compared with the classical 
collision-type heat-transfer coefficient. Substituting in (25) the density and 
temperature distributions from (27) yields 


« ~1.3X10-‘r ! /r e ! «c r* . ( 29 ) 

From (29) it can be shown that for a constant flux, the temperature gradient 
decreases as 


9T 

9r 


oc r 


-4 


B 


and may thus be responsible for the transition from rapid temperature increase 
[Serbu and Maier , 1966] to the isothermal temperature distribution of inter- 
planetary space far from the earth (see Figure 2 ). Because of these features 
of turbulent heat conduction, we adopt the turbulence process in the upper 
magnetosphere where the transition occurs from magnetosphere into inter- 
planetary space. 


Neglecting local and non-local heat sources at these altitudes, the energy- 
continuity equations for ions and electrons are 


V„ 




(i) 


+ V. 



= - L 


le » 


(30) 
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ALTITUDE (KM) 
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(31) 


• ( e ) . ( e ) 

v n 1 ii + V ± J ± V, =L 


ie, 


where the energy loss due to coulomb collisions between ions and electrons 
[ Brace, Spencer, Dalgarno , 1965] is given as 


Lie ~ 


(T i _ T e } 


(32) 


While the transverse turbulent heat conduction coefficients for electrons 

and ions are of the same order, the (reduced) parallel heat conduction 

coefficients differ by the factor (m./m )^ 2 . The downward electron heat flux 

i e 

parallel to the magnetic field into the lower atmosphere is thus the dominant 
process of cooling the upper part of the magnetosphere. The ions will lose 
their energy (which is conducted from interplanetary space into the mag- 
netosphere) mainly through coulomb collisions with electrons. Therefore, we 
expect the ion temperature to be larger than the electron temperature in this 
transition region. 


In order to estimate this temperature difference and the associated 
energy input from interplanetary space, we consider the term lie in (31) as 
a heat sink which must compensate the perpendicular ion heat transfer. This 
leads to 


J_ d_ 

r 2 3r 



- K 


N e 2 (Tj -T e ) 


T e 


w 


and after integration yields 


(33) 


15 


KN e 2 (Tj-Te) . 

T e ^ 


( 34 ) 





(i) 

J± 


(i) 


Jj-i 


3r 2 
r 3 - r, 3 


Equation 35 represents a relation between the ion-electron temperature 
difference, and the heat flux j ^ from the interplanetary space into the 

magnetosphere. Assuming that at a height of about r = 7r_ the heat flux 

(i) 1 

is almost attenuated (j =0), and adopting the values 


T e ~ 4 x 10 4 °K, 

K„ io" U 

f 

N ~ 10 cm -3 , 

we arrive with a reasonable temperature gradient of 0. 1 °K/km and a 
corresponding energy flux 

(i) st -io 2 

j x = - k bJT = - 10 erg/ cm sec 

at a temperature difference of 


Tj - T e ~ 10 4 # K. 

This is consistent with Serbu and Maier' s [1967] measurements. 

-10 i 2 

Evidently such heat fluxes in the order of 10 erg/cm sec are 

- 2~3 2 

extremely small compared with those energy fluxes of 10 - 10 erg/cm sec 
that arise from the non-local magnetospheric heating due to the escape of 
fast electrons from the ionosphere. 
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We conclude that the heat conduction between interplanetary space and mag- 
netosphere, though sufficient to explain the observed thermal non-equilibrium 
in the upper magnetosphere makes a negligible contribution to the energy 
balance of the magnetosphere. 

CONCLUSIONS 

The foregoing results make it possible to construct a model of mag- 
netospheric temperature distribution. Figure 2 shows a schematic equatorial 
temperature profile that extends from the bottom of the magnetosphere to 
interplanetary space. 

Three regions are defined: 

I. The lower magnetosphere, where the temperature distribution is 
determined by a normally high heat conduction due to high plasma densities 

3 4 -3 

{10 - 10 cm ), associated with heat received from the ionosphere. Here, 
the temperature profile is strongly affected by the observed latitudinal 
temperature increase up to 50 degrees geomagnetic latitude [ Brace, Reddy , 

Mayr , 1966] which is related to ion compositional features [Mayr, Brace, 
Dunham , 1967] in conjunction with the energy balance of the lower F region. 

II. The middle magnetosphere, where parallel heat conductivity can 

2 1 -3 

be largely decreased by low plasma densities (10 - 10 cm ), as previously 
discussed. The decrease in heat conduction and the consequent temperature 
increase along the field lines must be enough to counterbalance the temperature 
decrease with latitude observed at 1000 km above 50 degrees geomagnetic 
latitude [Brace, Reddy, Mayr, 1967] so that the temperature gradient increases 
with altitude in the equatorial plane. The energy input in this region may 
still be heating from the ionosphere. 
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III. The upper magnetosphere, where perpendicular heat transfer 
becomes significant. Here, turbulence and plasma instabilities f Bohm, 1949, 
Tsuda , 1967] enormously increase the perpendicular heat transport coefficient. 
The plasma components therefore can carry some heat energy from the mag- 
netopause and interplanetary space into the magnetosphere. Since the parallel 
heat conductivity is much greater for electrons than for ions, this incoming 
heat flux is transported mainly by electrons parallel to the magnetic field 
lines in the lower magnetosphere. The perpendicular heat conductivity of the 
ions is the same as that of the electrons; therefore, their heat energy is mainly 
transferred through coulomb collisions to the electrons. Thus ion temperature 

must be higher than electron temperature. With temperature gradients of 

— 1 — 10 2 
10 °K/km, energy fluxes of 10 erg/cm sec arise; these are large enough 

to explain the observed [ Serbu and Maier , 1967] temperature differences of 

4 

T^ - T e ~10 °K. This supports the concept of turbulent heat transfer in the 
upper magnetosphere. But it also proves that thermal coupling between 
interplanetary space and magnetosphere is negligible compared with the energy 
input from escaping fast ionospheric electrons which produce heat fluxes of 

- 2-3 .2 

10 - 10 erg/cm sec. 
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REFERENCES 

Bohm, D. , The characteristics of electric discharges in magnetic fields, 
edited by A. Guthrie and R. K. Wakerling, Chapter 2, Section 5, 
McGraw Hill, New York, 1949. 

Brace, L. H. , Reddy, and H. Mayr, Global behavior of the ionosphere at 
1000 km altitude, Geophys. Res. 72 . 265-183, 1967. 


18 


Brace, L. H. , N. W. Spencer, and A. Dalgarno, Detailed Behavior of the 
Midlatitude Ionosphere from the Explorer XVII Satellite Planet Space 
Sci. 13 , 647-666, 1965. 

Serbu, G. P. , and E. J. R. Maier, Thermal plasma measurements within 
the magnetosphere, Space Research VII, North Holland Publishing 
Company, Amsterdam, 1967. 

Spitzer, L. , Jr. , to Physics of fully ionized gases, Interscience Publishers, 
Inc. , New York, 1956. 

Tamao, T. , Temperature distribution in the magnetosphere, Report Ionosph. 
and Space Res. in Japan, 20, 1966, pp. 312-321. 

Taylor, Jr., H. A. , A. C. Brinton and L. R. Muenz, First results from 
OGO-C Ion Composition Experiment, Presented at the 47th Annual 
Meeting of AGU, March 1966. 

Carpenter, D. L. , and R. L. Smith, Whistler measurements of electron 
density in the magnetosphere, Rev. Geophys. 2 , 415-441, 1964. 

Geisler, J. E. , and S. A. Bowhill, Exchange of energy between the 

ionosphere and the protonosphere, Atm. Terr. Phys. 27 , 1119-1146, 
1965. 

Kaufman, A. N. , Dissipative effects, in: Plasma physics in theory and 

application, edited by W. B. Kunkel, McGraw Hill, New York, 1966, 
p. 92. 

Mayr, M. G. , L. H. Brace, and G. S. Dunham, On the ion composition and 
temperature in the topside ionosphere, J. Geophys. Res ., 1967. 


19 



Serbu, G. P. , and E. J. R. Maier, Low energy electrons measured on 
IMP 2, Geophys. Res. 71 , 3755-3766, 1966. 

Taylor, Jr., H. A., H. C. Brinton, and C. R. Smith, Positive ion 

Composition in the magnetosphere, obtained from the OGO-A Satellite, 
J. Geophys. Res., 70, 5769-5771, 1965. 

Tsuda, T. , Effective viscosity of a streaming collision free plasma in a 
weakly turbulent magnetic field, AGU-Meeting, Washington, D. C. , 
April 1967. 


20 


